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Ion distribution around a charged rod in one and two component solvents: 
Preferential solvation and first order ionization phase transition 
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Department of Physics, Kyoto University, Kyoto 606-8502, Japan 
(Dated: August 14, 2009) 

In one and two component polar solvents, we calculate the counterion distribution around an 
ionizable rod treating the degree of ionization a as an annealed variable dependent on its local 
environment. In the two component case, we take into account the preferential solvation of the 
charged particles and the short-range interaction between the rod and the solvent. It follows a 
composition-dependent mass action law. The composition becomes heterogeneous around a charged 
rod on a mesoscopic scale, strongly affecting the counterion distribution. We predict a first order 
phase transition of weak-to-strong ionization for hydrophobic chains. This transition line starts from 
a point on the solvent coexistence curve and ends at a critical point. The composition heterogeneity 
is long- ranged near the solvent critical point. 



I. INTRODUCTION 

Polyelectrolytes are extremely 

complicatedi£iM££I£&i i ii. Examples include bi- 
ological polyelectrolytes, like DNA, actin filaments or 
microtubules, and synthetic polyelectrolytes. In such 
charged polymers, the Coulomb repulsion among ionized 
monomers can induce a number of conformation changes 
of a chain. The Coulomb attraction among counterions 
and an ionized chain can result in condensation of coun- 
terions at large counterion contents (Manning-Oosawa 
counterion condensation). In practice, it is important 
that the phase behavior of polyelectrolytes strongly 
depends on the degree of ionization. For example, hy- 
drophobic polymer chains become soluble in water-like 
solvents with slight ionization. In this paper, we further 
investigate two complex aspects of polyelectrolytes, 
which have not yet been fully discussed. 

First, the ionization (or dissociation) process should be 
treated as a chemical reaction in many polyelectrolytes 
containing weak acidic monomers^ ' 12 ' 13 ' 14 ' 15 ' 16 ' 17 . The 
degree of ionization a is an annealed fluctuating variable 
governed by the mass action law and dependent on the lo- 
cal values of the counterion density and the composition 
(in mixture solvents). All these quantities depend on the 
electric potential self-consistently. Inhomogeneity of a 
appears on a chain 6 , but is crucial in structure formation 
and phase separation^. Even when a is nearly homo- 
geneous, it is a complex quantity dependent on various 
conditions. It can be small in solvents with low dielectric 
constant and can increase considerably in highly polar 
solvents. 

Second, complex effects are induced in polyelectrolytes 
when a second fluid component (cosolvent) is added to 
a water-like solvent. For example, precipitation of DNA 
has been widely observed with addition of (less polar) 
alcohol such as ethanol to wate r 18 i 19 ' 20 i 21 . Here the al- 
cohol added is excluded from condensed DNA 20,21 . On 
the contrary, with addition of zwitterionic species (which 
are more polarizable than water) , Flock et al^ observed 
a resolubilization of DNA in the presence of multivalent 



ions (as condensating agents). Baigl and Yoshikawa 23 
found that zwitterionic species increased stability of coil 
states of a single DNA molecule, which undergoes a dis- 
continuous phase transition between elongated coil and 
compact globule depending on the amounts of the sec- 
ond component and the multivalent ions (or cationic 
surfactant^) . 

To understand these mixture effects, relevance of the 
preferential solvation has been pointed out by experimen- 
tal group o 19 ! 20 ' 21 and by a theoretical group^. From 
our viewpoint, particularly important should be the ion- 
dipole interaction among charged particles and polar 
molecules 2 ^, which gives rise to the solvation (hydra- 
tion) shell composed of several solvent molecules (those 
of the more polar component in a mixture solvent) 
around each charged particle. The resultant solvation 
chemical potentials of ions typically much exceed the 
thermal energy ksT (per ion) and strongly depend on 
the composition in binary mixtures. It is decreased 
(increased) for hydrophilic (hydrophobic) ions with in- 
creasing the water composition in a mixture of wa- 
ter!- less polar component. Furthermore, the degree 
of ionization of polymers can strongly depend on the 
composition. Recently, including such solvation inter- 
actions, several theoretical groups have begun to in- 
vestigate the ion effects in electrolytes with mixture 
solvent o 25 i 27 ' 28 i 29 i 30 ' 31 , polyelectrolytes 1 -!, and ionic sur- 
factants at oil- water interfaces^ 2 -. 

In this work, we will demonstrate emergence of meso- 
scopically heterogeneous composition variations around 
a charged rod, which stem from the preferential solva- 
tion. If they are hydrophilic, the water-like component 
is enriched near a rod, even when the polymer backbone 
is hydrophobic. In such complex situations, the origi- 
nal concept of the Manning-Oosawa condensation is not 
available (or at least needs to be modified) to understand 
the counterion distribution. As a byproduct, we will pre- 
dict a first order phase transition between weakly and 
strongly ionized states of a rod. Here the ionization is 
assumed to be very weak for <pB = 0, but increase with 
increasing (ps- Our prediction is that the progress of ion- 
ization can occur as a discontinuous change. If a polymer 



2 



chain is in an expanded state in the weakly ionized phase, 
it should be more expanded in the strongly ionized phase. 
It is analogous to the prewetting phase transition of fluids 
on a planar boundary wal l 33 ' 34 i 35 . 

In Subsec.IIA, we will present a Ginzburg-Landau the- 
ory for one component solvents whose minimization gives 
a. In Subsec.IIB, we will analyze a and the counterion 
density n\ on the basis of some exact relations. In Sub- 
sec. IIIA, we will set up a Ginzburg-Landau model for 
two component solvents, which includes a, rii, and the 
composition <f> as fluctuating variables. In Subsec.IIIB, 
we will numerically examine these quantities in equilib- 
rium for various parameters without salt. In particular, 
we will treat a hydrophobic rod and hydrophilic counte- 
rions, where the dissociation sensitively depends on the 
ambient composition. In Appendix B, we will examine 
the ion distributions in one component solvent with salt. 
In Appendix C, we will present a simple thermodynamic 
theory of the solvation shell formation at small composi- 
tion 4>b of the more polar component . 



II. ROD IN ONE-COMPONENT SOLVENT 

We consider an ionizable polymer chain with a long 
persistence length in a one-component solvent. Its 
shape is a rod or an expanded coil. Necklace-like 
globules 3 ^^^ are outside the scope of this work, which 
appear in poor solvents with increasing ionization. As- 
suming low ion densities far from the rod, we neglect 
the formation of dipole pairs and ion clustering 4 and the 
free energy contribution from the charge density fluc- 
tuations (cx k 3 with k being the Debye-Hiickel wave 
number 

As a simple model 3 ^, a polymer chain is treated as a 
cylinder with radius b and length L 3> b The system is 
in a cylindrical cell with radius R^> b. We assume that 
all the quantities depend only on the distance r from the 
center of the rod neglecting the end effects. Here we ne- 
glect the discreteness of the charges along the chai n 11 ! 16 . 



A. Ginzburg-Landau free energy 

Let the density of ionizable groups on a rod be <jq and 
the degree of ionization be a in the range < a < 1. Wc 
assume homogeneity along the rod. Each ionized group 
has charge — e and the counterions are monovalent. The 
charge density along the rod is — ea with 



a$a. 



(2.1) 



per unit length. The mobile ions are distributed in the 
region b < r < R and < z < L. Their densities are 
written as rii(r) and their charges are Z^e {i = 1, 2, • • • ). 
In this work ni(r) denotes the density of the counterions 
from the rod plus the cations of the same species added 
as a strong salt (see Appendix B). Then Z\ = 1. The 



charge density in the region b < r < R is written as 

p = eJ2 z i n h (2-2) 

i 

where the summation is over all the mobile ions i = 
1, 2, • • • . We assume the overall charge neutrality, 



J drp(r) = Lea, 



(2.3) 



where J dr(- ■ • ) = 2ttL J b drr(- • • ) is the integral in the 
cell outside the rod. In the present one-component case 
the dielectric constant e is assumed to be a constant. The 
electric potential $(r) satisfies 



d 2 1 d \ ^ 

— r H $ = - 



dr 2 



dr J 



-47T/9. 



(2.4) 



We impose the boundary condition on the electric field 
E(r) = -d<&(r)/dr as 



E(R) = 0. 



(2.5) 



That is, there is no surface charge on the outer surface. 
See the review by Dobryinin and Rubinstein 7 for a the- 
ory in the case of nonvanishing E(R). With the aid of 
Eqs.(2.3) and (2.5), integration of Eq.(2.4) in the region 
b < r < R gives the boundary condition, 



E{b) = -2ea/eb. 



(2.6) 



Our theory remains invariant with respect to the shift of 
the potential $(r) — » $(r)+const. 

Hereafter we set the Boltzmann constant equal to 
unity. The Helmholtz free energy F of our system is 
divided in two parts as F = Fq + F^, where Fq con- 
sists of the entropic and electrostatic contributions and 
Fd is the free energy of dissociation. They are written 



as 



12,13,14,15,16,17 



Fa 
T 

Fa 
T 



J dr ni(\n(niV ) - 1) 



eE 2 
8t7t 



(2.7) 



= Lcr [a In a + (1 - a)ln(l - a)} + LA a. (2. 



In Fq, the volume vq is taken to be independent of i, 
which is allowable without loss of generality 3 ^. In Fd, 
Lao is the total number of the ionizable groups and Ao 
is the energy needed for ionization (see Eq.(2.14) for the 
dissociation constant in terms of Ao). The free ions can 
interact with polar segments on a chain, but such an 
interaction is neglected. 

In equilibrium we minimize the grand potential Q = 
F — J dr J^. \iiTii under the charge neutrality condition 
Eq.(2.3), where /ij (i = 1, 2, • • • ) are appropriate constant 
chemical potentials. Wc thus minimize 



Q = F 



dr - T\Zi)m - TXLa (2.9) 
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with respect to a and rii(r), where TA is the Lagrange 
multiplier. However, if we may replace A by A — fJ>x/T 
in F, we obtain Q = F without salt. 

The electrostatic energy F e = J dreE 2 /8ir changes 
with respect to small variations of a, rii, and e as 2 ® 



/E 2 
dr e$ ^2 z i^ n i - ~^~^ £ 



e$(b)6<r, (2.10) 



where use has been made of Eqs.(2.4)-(2.6). Here Se = 
for one component solvents, but Se will be nonvanishing 
for two component solvents. The minimum conditions 
dVl/da — Sfl/Srii = are rewritten as 



1 — Q 

n t (r) 



= exp 



e$(6) 



T 



n 0i exp 



-A-A 
/ e$(r) 



+ A 



(2.11) 
(2.12) 



where = Vq 1 exp(/ii/T) are constants. We notice 
that e$/T and A appear in the sum e$/T + A, so we are 
allowed to set A = by redefining e<f>/T + A as e$/T. 
Substitution of Eq.(2.12) into Eq.(2.4) yields the Poisson- 
Boltzmann equation. 

From Eq.(2.12) we have the relation ^ i [rin(n i w ) — 
fii + TXZi]rii = — p$. Further noting the relation 
f drp(r)[$(r) - $(6)] = 2F e and using Eq.(2.11) we find 
that U, assumes a negative minimum expressed as 

— = -a + a ln(l-a)-^-^, (2.13) 
TL v TL L y ' 

where N s = J dr^n rii — ah is the number of the ions 
added as a salt. In Eqs.(2.28) and (2.29) below, fl will 
be calculated explicitly without salt (N s = 0). 

As r — > b, we have n\{b) = n i exp[— e&(b)/T — A] for 
monovalent counterions. We thus obtain the mass action 
law (equation of ionization equilibrium), 



-ni(&) = K Q , 



1 

where Kq is the dissociation constant defined by 



Ko = n 01 e~ 



v n e 



IpW/T-Ao 



(2.14) 



(2.15) 



Since Eq.(2.14) yields a = 1/[1 + ni(b) / Kq], a increases 
up to unity for ni(b) -€i Ko. Note that Ko is a mea- 
surable quantity and should be invariant with respect to 
the choice of vo in Eq.(2.7)^. The mass action law of 
chemical reaction is often expressed in terms of the pH 
of the solution if the cations are protons^. In our case 
Eq.(2.14) holds for the counterion density ni(b) at the 
polymer backbone. 



B. Analysis in the salt- free case 

The Poisson-Boltzmann equation can be solved exactly 
without salt under the boundary conditions Eqs.(2.5) and 
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FIG. 1: Fi(r,a) = 27r£ s r 2 ni(r) vs ln(r/&) on a logarithmic 
scale for Ib<j = 0.8, 1, and 1.5 with M — 6 for one compo- 
nent solvent without salt. The slope of the curve at r/b — 1 
changes its sign at £b<j = 1. 



(2.6)^^1. The solution is parameterized by the normal- 
ized line density £b&, called the Manning parameter, and 
the radius ratio R/b. Here £b — e 2 /eT is the Bjerrum 
length. We examine how the degree of ionization a is 
determined for each given £bCo- ln this subsection we 
will give approximate results for large M . Some exact 
results are summarized in Appendix A. Furthermore, we 
will discuss the effect of salt in Appendix B. 



Counterion density and degree of ionization 



It is known that ni(r) cx r 2 and E(r 



roughly 



hold for r b. It is convenient to introduce the dimen- 
sionless function, 



Fi (r, cr ) = 2-Kl B r 2 m{r), 



(2.16) 



which behaves differently depending on whether a < a* 
or a > a* with a* being a critical line density, 

a* =£ B 1 /(1 + Ar 1 ). (2.17) 

Here we define the dimensionless parameter, 

M = \n(R/b), (2.18) 

which is assumed to be considerably larger than unity. 

In Fig,l, we show F\(r, a) for £b<J = 0.8, 1, and 1.5 at 
AI = 6 treating F\ as a function of u = ln(r/6). The area 
below each curve is the Manning parameter since 



r AI 

/ duFi(r, a) = £b<J. 
Jo 



(2.19) 



We notice that Fi (r, a) changes most drastically at r = b 
and most weakly at r = R. Its value at r = b is the 
normalized counterion density at the rod surface since 
F\(b,a) = 2ir£Bb 2 ni(b). To examine it, we introduce a 
scaling variable q by 



(£ B a - 1)M. 



(2.20) 



4 



i B <ro = 0.8 



3<T = 1.5 
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FIG. 2: Degree of ionization a (top) , normalized counterion 
density F\ (b, a) on the rod surface (middle) , and osmotic pres- 
sure II at r — R divided by Tn v = Ta/nR 2 (bottom) as func- 
tions of A in Eq.(2.22) for i B o = 0.8 (left) and 1.5 (right) in 
one component solvent without salt. 



At a = a* we have q - 
in Appendix A yield 



-M/(l+M) 



-1. The relations 



Fx(b,a) £* 2al B [alB - lfe 2q {-q » 1), 

= {ot B -\f (<?»1). (2.21) 

For |g| - 1 we have F x (b,a) ~ M~ 2 . Thus Fi(b,a) is 
very small for — q > 1 with increasing M, while it tends 
to be independent of M for q ^> 1. 

We now calculate the degree of ionization a from the 
mass action law Eq.(2.14). We introduce a normalized 
dissociation constant, 



A = 2ir£ B b 2 K . 
In terms of A, Eq.(2.14) is rewritten as 
aFi(b,a) = (1 - a) A. 



(2.22) 



(2.23) 



In Fig. 2, we plot numerical results of a, Fi(b,<r), and 
U/Tn p for £ B a = 0.8 (left) and 1.5 (right). Here we 




FIG. 3: Degree of ionization a vs A for Ibcq = 2, 3, and 
4 obtained from the asymptotic equation (2.26) in the limit 
R — > oo in one component solvent without salt. Here a — > 
1/^BO-Q as A — > 0, while a — > 1 for > 1. 



introduce the osmotic pressure at r = i?, 

n = Tni(i?) (2.24) 

where the Maxwell stress vanishes from Eq.(2.5). The 
n p = a/irR 2 is the counterion density for the uniform 
distribution. 

The first line of Eq. (2.21) gives 

(1 - £ B a a) 2 a 2 /(l - a) S A(R/b) 2 y2£ B a , 



(2.25) 



for — q = (1 — l B a)M ^> 1. The right hand side anoma- 
lously depends on R and is very small for R/b S> 1. On 
the other hand, the second line of Eq.(2.21) gives 



l B aoa - l) 2 a = A(l - a). 



(2.26) 



for q = (£ B (i — \ )M 3> 1. This cubic equation of a is the 
asymptotic equation independent of R, yielding a unique 
solution in the range (£ B (Jo)~ l < a < 1. Namely, a — > 
(^Bfo) 1 < 1 as A — > 0, while a — * 1 with increasing 
A. In Fig. 3, we show this limiting a obtained from 
Eq.(2.26) as a function of A for ^b<jo > h where the 
increase of a from to {£ b ctq)~ 1 in the narrow region 
< A < {b/R) 2 is not shown. 



2. Grand potential and effective polymer-solvent interaction 

In equilibrium without salt, the grand potential in 
Eq.(2.13) consists of three negative parts as 

Q/TL = -2a - r B x T e (£ B o) + a ln(l - a), (2.27) 

The electrostatic energy F e = J dreE 2 /8w is equal to 
TL(a + £g 1 J r e ). See Appendix A for the exact expression 
for the scaling function T e {s), which was derived by Naji 
and Netai^. (i) For small £ B a <§; M~ l the electrostatic 
part is negligible and il/TL = — a + <tq ln(l — a), (ii) For 
Af" 1 < £ B a «lwe have 



F e £3 TL£ B a M, 

Q/TL = -£ B a 2 M + <t ln(l 



a 



(2.28) 
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(iii) On the other hand, for i B a - 1 > M' 1 and £ B a < 
M, we find 



F e = TLl B l M, 

Sl/TL = -£ B X M + cr ln(l - a). 



(2.29) 



Here the electrostatic energy is nearly equal to that of a 
charged rod at a — £g without screening. This result 
suggests that a fraction of 1 — (isa)^ 1 of the counteri- 
ons are localized around the rod, as well as the osmotic 
pressure behavior in Eq.(A4). 

Polymer chains in water are often hydrophobic with- 
out ionization but can be hydrophilic with ionizations^. 
This means that the polymer-solvent interaction param- 
eter, written as x ps , is effectively decreased upon ioniza- 
tion. Its decrease A\ P s may be defined as follows. In 
the Flory-Huggins free energy density of polymer solu- 
tions, the interaction part is written as Tv$ Xps&p f° r 
small polymer volume fraction (f> p , where vq is the sol- 
vent volume. We map our polyelectrolytes system to a 
neutral polymer system. In the present case, the polymer 
volume is irb 2 L in the total volume V = irLR 2 so that 
ct> p = (b/R) 2 . We set 



Tv^&XvAp = si/v + il. 

Without salt we have II = Tn\(R) and 

" si 



A Xp 



irb 2 



TL 



KR 2 n l {R) 



(2.30) 



(2.31) 



Here in the brackets, the first term is given by Eq.(2.27) 
and the second term is equal to F\{R,(j)/21b from 
Eq.(2.16). We confirm the negativity of Ax ps as fol- 
lows, (i) For small £b& <C Af _1 we have Ax ps — 
{v /irb 2 )a a ln(l - a), (ii) For M^ 1 < l B a « 1 we 
have Axp S = (v /Trb 2 )[—£ B cr 2 + (r ln(l — a)], (iii) For 
i B a - 1 > M~ l and l B a < M, we have Ax ps = 
(v /irb 2 )Sl/TL, where SI is given by Eq.(2.29). 



III. ROD IN TWO COMPONENT SOLVENT 

We next consider a charged rod in a nearly incom- 
pressible two component solvent such as a mixture of 
water+alcohol or water + organic solvent in the same ge- 
ometry as in the one component case. The volume frac- 
tion of the water-like component is <p(r) and that of the 
less polar component is 1 — </>(r). The molecular volumes 
of the two components (the inverse densities of the pure 
components) are assumed to be given by a common vol- 
ume vq = a 3 , which may be equated with vq in Eq.(2.7). 
Hence <j> will also be called the composition. We assume 
small ions and neglect their volume fraction. As in the 
previous section, n\(r) is the density of the counterions 
plus the cations of the same species added as a salt. 

We suppose a polymer chain with a hydrophobic back- 
bone (like polystyrene) and ionizable groups attached to 
itl. With addition of water, we first need to consider 



the formation of solvation (or hydration) shells composed 
of several water molecules around hydrophilic charged 
particles 26 . This can occur even at very small water 
content^, as will be discussed in Appendix C. 



A. Ginzburg-Landau theory 

For not very small bulk water composition <p Bl the 
water composition <fi(r) varies on mesoscopic scales. To 
describe such situations, we present a Ginzburg-Landau 
theory with the gradient free energy of the composition 36 . 



1. Free energy including solvation interaction 

The total free energy is composed of three parts as 
F = F + F d + AF, where F is given in Eq.(2.7) and 
F d in Eq.(2.8). Notice that A in Eq.(2.7) should be 
replaced by A in Eq.(C7) to account for the solvation 
shell formation. Hereafter we will redefine Ao as the right 
hand side of Eq. (C7) . Then TAo in the following is the 
renormalized dissociation energy per counterion. 

Assuming the homogeneity along the rod, we write the 
additional contribution A_F in the form, 



AF 



dr 



k{<t>) 



T 

2TibL^f(j) s 



+ f|V0| 2 
LAidoa^, 



(3.1) 



The space integral is in the cell, b < r < R and < z < L. 
The free energy density /o(</ > ) is taken to be the Bragg- 
Williams form, 



T 

/o = - 



> In <t> + (1 - 4>) ln(l - 4>) + X<t>{l ~ <t>)l (3.2) 



where x is the interaction parameter dependent on T 
and its mean-field critical value is 2 in the absence of 
ions. The coefficient C of the gradient part is a positive 
constant of order a" 1 and will be taken to be 3a _1 in our 
numerical analysis. The coupling terms (oc gi) arise from 
the ion-dipole interactions among the ions and the polar 
solvent molecules. In the second line of Eq.(3.1), the term 
proportional to 7 arises from the short-range interac- 
tion between the rod surface and the solvent molecules 3 - 3 -,, 
while the term proportional to Ai represents the solva- 
tion interaction between the surface charge and the sol- 
vent molecules. We neglect the interaction between the 
mobile ions and the uncharged monomers (which can be 
important for polar monomers^). Hereafter, 



(3.3) 



is the surface value of the composition (outside the sol- 
vation shells under the condition (C8)). These molecular 
interactions of the rod are characterized by the two con- 
stants 7 and Ai, which can be either negative or positive. 



6 



In Fq in Eq.(2.7), the dielectric constant e is assumed 
to depend on the composition in the linear form 4 ^, 



e{4>) = s + ei 



(3.4) 



where eo + Si is the dielectric constant of the first compo- 
nent and £o is that of the second component. Thus e{<j>) 
is inhomogeneous. The electric potential <j>(r) satisfies 



--^-re(</>)-^-$ = -4tt P , 
r ar ar 



(3.5) 



under the boundary conditions E(R) = and 

E{b) = -2ea/e((j) s )b. (3.6) 

2. Equilibrium and composition- dependent mass action law 

We minimize the grand potential Q defined as in 
Eq.(2.9) with F = F + F d + AF. The counterparts 
of Eqs.(2.11) and (2.12) read 



a 



= cxp 



1 - a 
m(r) = n 0i exp 



A + 



e$(6) 
z e<&(r) 



A + Ai 



, (3.7) 
(3.8) 



where noi = v ~ 1 e tJ,i ^ T ~ ZiX . As in Sec.l, we may set A = 
without loss of generality. From these equations the mass 
action law follows as 



1 



-m(b) = K{ct >s ). 



(3.9) 



The dissociation constant K{<p s ) depends on the surface 
composition <f> s as 



K(<j> s 



n xe 



-A + (Ai+gi) 



(3.10) 



With increasing <j> 8 , K{<j> s ) increases (decreases) for posi- 
tive (negative) Ai + <?i. We are not aware of experimental 
data on the composition dependence of a for polymers 
in mixture solvents. On the other hand, strong compo- 
sition dependence has been reported in dissociation of 
weak acids in various aqueous mixture o 42 ' 43 . 

The functional derivative h — 5F/5<fi at fixed ni and a, 
is the chemical potential difference of the two components 
divided by t>o and is homogeneous in equilibrium. From 
Eqs.(2.10) and (3.1) we obtain 

h = f Q {J>) - TCV 2 ^ - TYa ini - ^E 2 , (3.11) 



where 4>b — 4>{R) is the value oi <j) dX r — R. As in 
Eq.(2.13) some calculations give 



n 

TL 



2ir / drr 

Jb 



-a + cr m(l - a) - — - — . 



2irlrf(p s 



In the first line we write 4>' = d<p/dr and define 

/o(0)=/oW-/o(0B)-ftW-^fl). 



(3.13) 



(3.14) 



The first line of Eq.(3.13) gives the compositional contri- 
bution. The second line is of the same form as the right 
hand side of Eq.(2.14) with N s being the number of ions 
added as salt. 

Small variations of (j>, a, and n% yield an incremen- 
tal change of £1, written as SQ. With the equilibrium 
conditions (3.5), (3.7), (3.8), and (3.11), its linear terms 
proportional to 5a and Srii vanish. There remain surface 
parts of the form, 



H = 27rC[(n^) r=a 



{r<t>'54>)r=b] 
-(2tt67 + Aia)8(j>(b). 



(3.15) 



We set 5il — for any boundary variations 8<f>(b) and 
5cj)(R) to obtain 

<j>'(R) = 0, (3.16) 
C<f>'(b) = -7 - (Ai/2tt6)ct, (3.17) 

which are the boundary conditions on <p i n solving the 
equation h =const. in Eq.(3.11). In our theory the con- 
dition at r = b in Eq.(3.17) depends on a. For positive 
(negative) 7 + Aicr/27r6, the first component tends to be 
attracted to (repelled from) the rod. A similar bound- 
ary condition has been used in the gradient theory of the 
wetting transition on a planar wall^ 3 .. 

We use the gradient free energy. It is usually derived 
from the gradient expansion of the two body van der 
Waals interactions and is justified near the critical point. 
In the following, however, we will present numerical re- 
sults in the composition range 0.15 < 4>b < 0.53, impos- 
ing the boundary condition (3.17) in the gradient theory, 
where 4>b is the composition far from the rod. When the 
rod radius b is of order a and <p changes steeply around 
the rod, a density functional theory without the gradient 
expansion should yield more reliable result o 34 i 35 . 



8tt 



where f' {<j>) — dfo((/))/d(f> and E = —d^/dr. In equilib- 
rium ft, is a homogeneous constant. 

With these results it is convenient to redefine the grand 
potential for two component solvent as 



fl = F 



dr 



^ mni + /o(0s) + h(d> - (f) B ) 



, (3.12) 



3. Solvation interaction 

In the original Born theory^ 4 -, the solvation chemical 
potential was given by n % so \{4>) = Zfe 2 /2e(<j))Ri for ion 
species i, where Zie is the charge, e(</>) is the solvent di- 
electric constant dependent on <fi, and Ri is a microscopic 
length called the Born radius. However, this formula is a 
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crude approximation and is not applicable to hydropho- 
bic ions. In this paper we assume the form 

mLiW = ^ o1 (0)-^. (3.18) 

The solvation contribution to the free energy density is 
given by J2i f x loi( ( t > ) n i^ yielding the terms proportional 
to gi on the right hand side of Eq.(3.1). For a mixture 
of water and oil, for example, Tgi is the difference of 
the solvation chemical potential in oil and that in water. 
Therefore, > for hydrophilic ions and gi < for 
hydrophobic ions in aqueous mixture solvents. 

The difference of //* ol between the coexisting phases co- 
incides with the Gibbs transfer energy AG£ r (per ion) in 
electrochemistr y 39 ! 45 . It becomes Tgi in strong segrega- 
tion in our theory. Data of AG£ r for water+nitrobenzene 
at T ~ 300K suggest gi ~ 15 for monovalent ions such 
as Na + and gi ~ — 15 for tetrarphenylborate BPhJ . For 
water+alcohol, should not be smaller, since the di- 
electric constant of nitrobenzene is 35 and that of alcohol 
is smaller (24 for ethanol). The solvation coupling is thus 
very strong in aqueous mixtures. 

4- Analysts of behavior of composition deviation 

The composition <p(r) tends to a constant <pB far from 
the rod. We treat 4>b as a control parameter. With a salt 
added, this behavior is obvious if the Debye length k~j^ 
is shorter than the system radius R and longer than the 
correlation length £ of the composition. In such a case 
Eq.(3.11) gives <f>(r) — <\>b cx e~ kDT far from the rod. 

Without salt, however, the decay is algebraic. To show 
it, we linearize Eq.(3.11) with respect to the composition 
deviation 6<f> — <f> — 4>b to obtain 

(l-£ 2 V 2 )^ Xs .9im + §^£ 2 , (3.19) 

where we assume \S<j>\ <C 4>b at any r. We define the 
correlation length £ and the susceptibility \ s as 

£ = [TC/fZ(<p B )} 1/2 , (3.20) 
Xs = T/fH(<t> B )=e/C, (3.21) 

where fg{4>) = T[l/0(1 - <p) - 2x]/v from Eq.(3.11). In 
this paper, we set C = 3/ a; then, Xs ~ a£ 2 - The length 
£ remains of order a away from the criticality, while it 
grows near the criticality. 

The two terms on the right hand side of Eq.(3.19) 
roughly behave as r~ 2 from the results for one component 
solvent. That is, if Ft = 2n£ B0 r 2 n 1 (r) in Eq.(2.16) (or 
in Eq.(3.29) below) is of order unity, the two terms are 
of order Xs9i/2Tr£Bb r2 and Xs£i/2nEo£Bb'r 2 , respectively. 
We define two Bjerrum lengths as 

(so = e 2 /e T, £ Bb = e 2 /e B T, (3.22) 

at e = £o an d e = £(4>b) = £b, respectively. Thus Isb = 
^bo/(1 + ^i < / , b/£o)- Furthermore, we are interested in the 



strong solvation case gi ^> 1 with Si/eq ~ 1, where the 
first term dominates over the second in the right hand 
side of Eq.(3.19). Then 8<f)(r) decays for r — b > £ as 

54>(r) ~ { gi at 2 /2nl Bb )r- 2 . (3.23) 

Here we are assuming \5<fi\ < <j)B, which holds even for 
r-~£if 5l a/2^B6<l. 

5. Effective charge density and modified Manning law 

For large r far from the rod, 4>{r) tends to a constant 
(f)B and ni(r) obeys the nonlinear Poisson-Boltzmann 
equation even in the salt-free case. We introduce a char- 
acteristic length r c in the range b < r c -C R. For r larger 
than r c , we may set e((f>) = £b with sb — £(</ ) s) and the 
potential <&(r) obeys 

-e B V 2 $ ^ 47ren exp(-e$/T), (3.24) 

where we set Z\ = \ and no — noie -A . We compare $(r) 
and n\{r) in our system approximately obeying Eq.(3.24) 
far from the rod and those obtained as the solution of the 
nonlinear Poisson-Boltzmann equation, 

- £bV 2 $pb = Anen\ B , (3.25) 

where n\ B (r) = no exp(— e$ps (r)/T) is the correspond- 
ing counterion density. We set <j> PB (R) — at r = R and 
impose the boundary condition at r = r c as 

*ra(r<0 - *'(r c ), (3.26) 

where $ PB (r) = d^ PB {r)/dr and $'(r) = d<&(r)/dr. 
From Eq.(3.25) $(r) = $pbM and n\(r) = n^ B (r) far 
from the rod r > r c , but significant differences can arise 
in the region r < r c . 

In terms of $pb we define the effective charge density 
a e g on the rod by 

tfeff - e B b<P' PB (b)/2e. (3.27) 

From Eq.(2.6) the real charge density a is given by a = 
e{<j> a )W{b)f2e. From Eqs.(3.5), (3.24), and (3.25) the 
apparent excess charge density Ac = a — cr e g on the rod 
is expressed as 

Aa = 2nn f° drr[e" e *M/ T - e" e * PB «/ T ] 

Jb 

= ^- q £ du[F x {u)-Fr{u)]. (3.28) 

In the second line, we set u — ln(r/fr) with u c — ln(r c /6) 
being the upper bound and define 

F? B = 2nr 2 e B0 n^ B , (3.29) 

as in Eq.(2.16). These quantities are functions of u. In 
the bottom plates of Fig. 4, we shall see that u c may be 
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pushed to infinity since the integrand in the second line of 
Eq.(3.28) vanishes for large u. Here the effective Manning 
parameter is <J c tt(-Bb, where £ B b is defined in Eq.(3.22). 
In terms of n° ff = a e g/TrR 2 the asymptotic law for the 
osmotic pressure II in Eq. (A4) is changed to 



n 



Tn cS 



(1 - £ B b<T C ti/2) (^B&CTcff < 1) 

{£ Bb a eS > 1). (3.30) 



Tnf/2£ Bb a eS 



6. Strong deformations of the counterion distribution 

We examine the conditions of strong attraction of the 
counterions due to a composition change around the rod. 
We assume that <f> B is not very small and £ is not much 
separated from a for simplicity. If gi 3> 1, 8(f) is of order 
g\aj2-K£ Bb at r ~ £ from Eq.(3.23). Since n\ oc e 91 * from 
Eq.(3.8), an appreciable attraction of the counterions is 
induced due to the composition change when 



gla/2w£ B b > 1. 



(3.31) 



condition can be realized 
b). 



For large g% ^> 1 the above 
while 8(f> remains small (<C 

Also </)'(b) satisfies the boundary condition (3.17). This 
yields a contribution to 5(f>(r), written as (8(f)) b (r). It 
obeys (1— £ 2 V 2 )(60)& = in the bulk, so it decays rapidly 
far from the rod. If £ < &, it is approximately written as 



(ty)b(r 



7' 



Aicr 
"2^6 



^exp 



(3.32) 



If the above composition change is positive, the counte- 
rions are significantly attracted to the rod when 



. gi (7 + A 1 a/2 7 r6)£/C , > 1. 



(3.33) 



If it is negative and its absolute value exceeds unity and if 
the condition (3.31) does not hold, the counterions should 
be repelled from the rod. 

The criterions (3.31) and (3.33) are crude ones based 
on many assumptions. In particular, F\{b) ~ 1 has been 
assumed in deriving Eq.(3.31) and the unknown a = o^cx 
is contained in Eq.(3.33). Setting up general criterions 
is at present difficult, because many of the parameters 
strongly affect the counterion distribution. Neverthe- 
less, we recognize that the counterion distribution can 
be changed dramatically even for a slight change of the 
composition in the strong solvation condition |<?i| 3> 1. 



B. Numerical results for two component solvent 
without salt 



We present some numerical results for a water-oil sol- 
vent without salt. In all the examples to follow, we nu- 
merically solve Eqs.(3.5), (3.7), (3.8), and (3.11) under 



given boundary conditions by setting 

b = 2a, R/b = e 6 = 403, 

cr = 7r/2a, C = 3/a, Ai = 7, 

ex = 2e , £bo = 12a/ir = 6/<j , 

where £bq is defined in Eq.(3.22). The solvation param- 
eter gi will be either of 5, 7, or 10. Note that g\ can 
be larger in real situations, as discussed below Eq.(3.18). 
The parameter A = 2ir£ Bob 2 K in Eq.(2.22) will be cho- 
sen to be small. As a result, little ionization occurs with- 
out composition variations near the rod. In our simula- 
tion the counterion density far from the rod is so small 
such that the solvent phase diagram is unchanged 2 - 8 . 



Attraction of water and counterions to a hydrophobic rod 




1 1.5 

ln(r/6) 



FIG. 4: Numerical results for a mixture solvent in the salt-free 
case for \ = 1.4, <f>B = 1/3, gi = 5, and A = 0.0096 (on the 
point (+) in Fig. 6). Top plates: composition cf>(r) (left) and 
normalized counterion density voni(r) (right) for a 2 7 = —0.6 
and -0.2. Bottom plates: Fi(u) and F[ B (u) in Eq.(3.29) 
vsu = ln(r/&) for a 2 7 = -0.6 (right) and -0.2 (left). The 
area between these two functions (shaded) is equal to the 
normalized excess charge density ^boAct from Eq.(3.28). 



The interaction between the solvent and the rod in 
Eq.(3.1) leads to the boundary condition (3.17) in our 
gradient theory. For 7 < (7 > 0), the rod repels 
(attracts) the water component and is hydrophobic (hy- 
drophilic) without ionization. Even for 7 < 0, how- 
ever, the right hand side of Eq.(3.17) changes from pos- 
itive to negative with increasing the degree of ioniza- 
tion a and the condition (3.33) can eventually hold for 
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7I < Ai<r/27r6. Then the rod becomes effectively hy- 
drophilic. To examine the resultant preferential adsorp- 
tion of the water-like component, we introduce 



T = 2tt 



drr[(j)(r) — <f>i 



(3.34) 



where <pg is the value of <f>(r) at r = R. 

Figure 4 illustrates such a changeover in the one phase 
region for (a) a 2 7 = —0.6 and (b) a 2 j = —0.2. We set 
X = 1.4, 4> B = 1/3, gi = 5, and v Q K Q = 10~ 4 . Then 
A = 2Tr£ B ob 2 K = 0.0096 from Eq.(2.22), £ = 1.32a 
from Eq.(3.20), and v K = 10~ 4 e 12 ^ from Eq.(3.10). 
In the upper left panel of Fig. 4, we present the com- 
position profile <f>{r) near the rod, which is repelled 
for (a) and attracted for (b). However, as shown in 



Eq.(3.23), 



has a positive tail 



■) with- 



out salt, giving rise to a positive logarithmic contribu- 
tion - (gia£, 2 /£ B0 ) ln(i?/0 to the integral T in Eq.(3.34). 
Due to this singular contribution we obtain T = 0.706a 2 
for (a), while T = 3.99a 2 for (b). In the upper right 
panel, the difference of ni(r) between the two cases is 
large and n\{r) is written on a semi- logarithmic scale. 
Here, a = 0.369, <j)(b) — 4>b = —0.101, and wo«i(o) = 
2.73 x 10~ 3 for (a), while a = 0.621, <j>{b) - <j> B = 0.214, 
v ni{b) = 4.19 x 10~ 2 for (b). The lower panels of Fig. 4 
display the two functions F\{u) and Ff B {u) in Eq.(3.29) 
for the two cases. The excess charge density in Eq.(3.28) 
is negative as Act = — 0.02/^b = — 0.03/^_bo f° r ( a ), but 
is positive as Act = 0. 33 /£sb = 0.56/^so for (b). 

We then check the criterions (3.31) and (3.33). The left 
hand side of Eq.(3.31) is 1.73, while that of Eq.(3.33) is 
-0.609 for (a) and 0.755 for (b). However, F x {b) = 0.109 
for (a) and 1.54 for (b), so the criterion (3.31) is not 
justified for (a) (see the discussion above Eq.(3.22)). 



2. First order phase transition of ionization without salt 
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FIG. 5: Right: Phase diagram in the <j>B~X plane, where gi = 
10, Kovo — 8 x 10~ 6 , and a 2 j = —0.4 for water-oil solvent 
without salt. Left: Normalized grand potential aQ./LT on 
lines (A)-(E). For (A)-(C) there are two branches of weak and 
strong ionization and the equilibrium is given by the lower 
branch, while for (D) and (E) there is only one branch. 
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FIG. 6: Degree of ionization a vs 4>b without salt along the 
lines (A)~(E) in Fig.5, where gi = 10, K v = 8 x 10" 6 , 
and o 2 7 = —0.4. The discontinuity (broken line segments) 
vanishes at the critical point of ionization. 



With increasing 4>b we predict a first order phase 
transition from weak to strong ionization. We sup- 
pose that the backbone is hydrophobic, the ionization 
is weak in the pure second component, and the counteri- 
ons are strongly hydrophilic. Hence, in Figs. 5-9, we set 
a 2 7 = -0.4 x 10~ 6 , vqK q = 8, and g x = 10. 

We vary cf>B at fixed \ on the lines (A)-(E) in the left 
panel of Fig.5, where x = 2.2 (A), 1.7 (B), 1.3 (C), 1.16 
(D), and 1 (E). For x > X™? there are two branches of 
weak and strong ionization around a first order transition 
line expressed as 

4>b = <^ a (x)- (3.35) 

It starts from a point given by {4>b, X) — (0.1736, 2.39) on 
the solvent coexistence curve and ends at a critical point 
given by (<f> B ,x) = «,X C ") = (0.415,1.16), where 
= < / ) B a (x crl )- l n t ne right panel of Fig.5, the grand 
potential Q in Eq.(3.13) is calculated, which is lower on 
the equilibrium branch than on the metastable one. In 



Fig. 6, we show a vs 4>b on the lines (A)-(E), where (A)- 
(C) pass through the transition line, (D) passes through 
the critical point, and (E) is a supercritical path. This 
ionization transition is analogous to the prewetting phase 
transition on a wal l 33 ' 34 i 35 , where a first order phase tran- 
sition line also starts from the coexistence curve ending 
at a critical point. 

In the left panels of Fig. 7, we show the charge densi- 
ties ct and CT c ff multiplied by tsb on the lines (A) and (C), 
where Isb is defined in Eq.(3.22) and a e g by Eq.(3.27). 
In these cases, £_BbC increases but iBb&eS decreases at the 
transition from weak to strong ionization, The counteri- 
ons are more strongly attracted to the rod in the strongly 
ionized state than in the weakly ionized state (see Fig. 8 
below). The right panels of Fig. 7 display the counterion 
density n\{R) = II/T on the outer surface as a function 
of 4>b- Counter- intuitively, n decreases discontinuously 
at the transition with increasing a. The jump of n at 
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FIG. 7: Charge density a and effective charge density a e g 
in Eq.(3.19) on the rod around the transition in units of t^ b 
(left). Counterion density ni(R) = Tl/T on the outer surface 
in units of Vq 1 (right). These are plotted as functions of <f>B on 
the line (A) (upper plates) and on the line (C) (lower plates) 
in Fig. 5. In the weakly ionized phase, Isb^ is smaller than 
unity for (A) and larger than unity for (C), leading to larger 
discontinuities for (A) than for (C). 
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FIG. 8: Profiles without salt at the two points (al) (weakly 
ionized state) and (a2) (strongly ionized state) on the line 
(A) in Fig. 5. Top plates: 4>(r) (left) and wofii(r) (right) vs 
r/b. Bottom plates: Fi and F[ B vs ln(r/6) for (al)(left) and 
(a2) (right). From Eq.(3.28) the area of the gray region is 
equal to ~£boA<j for (al) and to IboAct for (a2). 



the transition is of order 8% on the line (A) and 1% on 
the line (C) in accord with the modified Manning lim- 
iting law (3.30). Here the effective Manning parameter 
^BbCcS is smaller than unity on the line (A) and larger 
than unity on the line (C). 

In Fig. 8, we show the spatial profiles of </>(r), ni(r), 
Fi(r), and F^ B (r) at the two points (al) and (a2) on 
the line (A) in Fig. 5. The left hand sides of Eqs.(3.31) 
and (3.33) are equal to 5.83 and —0.877 for (al) and to 
6.13 and 2.53 for (a2), respectively. Here, a = 0.221, 
(f)b = 0.143, v ni(b) = 3.20 x 10" 4 , and T = 0.706a 2 at 
(al), while a = 0.998, cj) b = 0.982, u ni(6) = 0.245, and 
T = 3.99a 2 at (a2). The upper left panel of Fig. 8 shows 
that the water component is considerably depleted from 
the rod at (al) but it covers the rod almost completely 
at (a2). The upper right panel of Fig. 8 indicates that the 
counterions are more accumulated around the rod at (a2) 
than at (al) such that they are more depleted far from 
the rod at (a2) than at (al). In the lower panels of Fig. 8 
the excess charge density Ac in Eq.(3.28) is — 0.003/f_Bb 
at (al) and 3.32/£ Bb at (a2). 
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FIG. 9: Change of the effective polymer-solvent interaction 
parameter Ax P s due to ionization for a two component solvent 
without salt along the paths (A)-(E) in Fig. 5. 

In Fig. 9, we plot the change of the effective polymer- 
solvent interaction parameter A\ ps in Eq.(2.31) along 
the paths (A)-(E). Its negativity indicates that the sol- 
vent quality becomes effectively better with ionization 
and adsorption. In the brackets of Eq.(2.31), the second 
term ttR 2 ii{R) is at most 10% of the first term fl/TL 
and its discontinuity at the transition is very small. 

We comment on the grand potential fl in Eq.(3.13). 
At the transition, the composition part (the first line) in- 
creases due to the layer formation, while the dissociation 
part (the first two terms in the second line) decreases. 
The change of the electrostatic part is much smaller than 
these changes in the present case. At the transition point 
on the line A, these three parts (multiplied by a) are given 
by (0.848, -0.775, -1.36) and (11.0, -11.0, -1.28) in the 
weakly and strongly ionized states, respectively. 

The first order phase transition can occur over a wide 
range of the parameters both without and with salt. In 
Fig. 10, the first order phase transition lines are displayed 
for various g\ and 7 without salt, which start from the 
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FIG. 10: First-order phase transition lines of ionization in the 
4>b~X plane for various parameter values without salt. Solid 
lines: a 2r y = —0.4, —0.6, and —0.8 with gi = 7 and Kovo = 
8 x 10~ 6 . Broken lines: a 2 ^ = -0.2,-0.4, and -0.6 with 
gi — 5 and Kovo = 10 -4 . Each line starts from a point on 
the solvent coexistence curve and ends at a critical point with 
decreasing \ m the bulk one-phase region. 



solvent coexistence curve to end at an ionization critical 
point. These lines are markedly enlarged with increasing 
the hydrophilic solvation strength (<?i > 0) and/or the 
rod hydrophobicity (7 < 0). The transition with salt 
will be examined in future. 



3. Profiles near the solvent coexistence curve without salt 

In the following we examine the profiles of <j)(r) and 
ni(r) close to the water-poor branch of the solvent co- 
existence curve (4>b < 1/2)- We vary (\>b and x fixing 
the other parameters as in Figs. 5-9. Namely, g\ = 10, 
a 2 7 = -0.4, and v K = 8 x 10~ 6 . 

Figure 11 presents the profiles at two points, {4>b, X) = 
(0.1730,2.392) and (0.1740,2.389), between the transi- 
tion point (4>b,x) — (0.1736,2.390) on the coexistence 
curve. We recognize marked jumps at the transition. 
That is, at the weakly ionized state at 4>b — 0.173, 
we have a = 0.202, cr = 1.2W"J = 0.900^, and 
r = 0.565a 2 , while at the strongly ionized state at 



4> B = 0.174, we have a = 0.999, a = 5.99^ = 4.44^, 
and r = 27.2a 2 . The first component remains repelled 
around the rod in the weakly ionized phase but is much 
attracted around it in the strongly ionized phase. 

In the upper panels of Fig. 12, we set 4>b = 0.45 to 
obtain T/a 2 = 86.6, 178, and 497 for X = 1-9,2.0, and 
2.0067, respectively. The point {<f) B , x) = (0.45,2.0067) 
is on the coexistence curve. In the lower panels, we set 
4> B = 0.5 to obtain T/a 2 = 88.9,211.0, and 8691 for 
X = 1-97, 1.990, and 2.0, respectively. A shoulder in </>(r) 
at (4>b,x) — (0.45,2.03) (left in the top) represents a 
water-rich layer varying gradually. At the solvent critical 
point, the excess deviation 8(f>(r) — 4>{r) — cf>B decays as 
5<f>(r) ~ r _1 for r/a < 10 and 54>(r) ~ r~ 2 for 10 < 
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FIG. 11: Composition (f>(r) and normalized counterion density 
von\(r) on the coexistence curve for {<j)B,x) — (0.1730, 2.392) 
in the weakly ionized phase (bold lines) and (0.1740,2.389) 
in the strongly ionized phase (dotted lines). The line of the 
first order ionization phase transition starts between these two 
points as in the right panel of Fig. 5. 
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FIG. 12: Composition 4>(r) and normalized counterion density 
voni(r) near the solvent coexistence curve for 0.45 (top), and 
0.5 (bottom). Values of \ and £ are given within each panel. 



r/a < 50. In all these cases, the dissociation is nearly 
complete or a = 1. 

Marcus et obtained a water-rich layer separated 
from the water-poor bulk region by a sharp boundary 
around a rod or a sphere. In our scheme, a layer with 
an interface follows in the low density limit of the coun- 
tcrions (realized for very large R) and under the condi- 
tions 7 = Ai =0 (which yields the boundary condition 
<f/(b) = from Eq.(3.17)). In addition, if the bulk re- 
gion is metastable (inside the solvent coexistence curve), 
a layer around a chain can trigger phase separation. This 
is analogous to the nucleation process around a charged 
particle in metastable polar fluids^. 
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IV. SUMMARY AND REMARKS 

In this paper, the ion distributions have been examined 
around a charged rod when the degree of ionization a is 
a fluctuating quantity. In mixture solvents, the effect of 
the preferential solvation has been investigated. 

In Sec. II, a long ionizable rod in a one component 
solvent has been treated, where the dissociation process 
gives rise to the free energy contribution Fd in Eq.(2.8). 
Minimization of the grand potential fl in Eq.(2.9) has 
then yielded a obeying the mass action law and the 
charge distribution n\{r) obeying the Poisson-Boltzmann 
equation. In Subsec.IIB, we have examined the counte- 
rion density at the rod surface without salt on the basis of 
the exact results 3 -!. For considerably large M = \n(R/b), 
a is determined by Eq.(2.25) for (1 — oIb)M » 1 and by 
Eq.(2.26) for (ol B ~ l)M > 1. All the quantities sensi- 
tively depend on the parameter ^4(oc Kq) in Eq.(2.22) 
as in Fig. 2. In the limit M 3> 1, a tends to unity 
in the former and to a well-defined limit in the range 
(^bCo)" 1 < a < 1 in the latter as in Fig. 3. 

In Subsec.IIIA, we have generalized our theory in Sec. II 
to describe the ionization equilibrium in mixture solvents. 
The additional free energy is AF in Eq.(3.1) for the com- 
position </>, which includes the solvation couplings with 
ions (oc gi) and the ionized monomers (oc Ai). The 
Manning limiting law for the osmotic pressure (A4) is 
modified to Eq. (3.30). Though still fragmentary, Sub- 
sec. IIIB has presented numerical results without salt, 
where the solvent consists of a water-like component and 
a less polar component. The counterions and the charged 
monomers are hydrophilic with g\ > and Ai > 0. 
Even if a rod is hydrophobic with 7 < 0, it becomes 
effectively hydrophilic with ionization as in Fig. 4. We 
have found a first order phase transition of ionization for 
hydrophilic counterions around a hydrophobic rod as in 
Figs. 5-8. At the transition from weakly to strongly ion- 
ized states, the number of the counterions increases, but 
the osmotic pressure has decreased in our examples as 
in Fig. 7. The polymer-solvent interaction parameter Xps 
decreases upon ionization as in Fig. 9. We have examined 
the composition and counterion profiles at the crosspoint 
of the ionization transition line and the coexistence curve 
in Fig. 11 and near the coexistence curve in the strongly 
ionized phase in Fig. 12. The adsorption T of the compo- 
sition in Eq.(3.34) much increases near the coexistence 
curve in the strongly ionized phase. The adsorption is 
long-ranged near the solvent criticality. In Appendix B, 
we will add a salt, where a is a decreasing function of the 
salt density. In Appendix C, we will examine the solva- 
tion shell formation at small content of a polar solute. 

We make further remarks on the first order phase tran- 
sition of ionization of an ionizable rod. Here we suppose 
that a polymer chain can be in an expanded coil state in 
a weakly ionized state and in a more expanded state after 
the transition without coil-globule transition. We have 
calculated inhomogeneities perpendicular to the chain, 
but those along the chain should also be important and 



can well alter the nature of the transition. 

Finally, we mention previous and proposed experi- 
ments, (i) Many authors have studied the conformation 
of a neutral polymer near the solvent critical poin t 48 i 49 i 50 . 
The effect of the critical fluctuations is much more en- 
hanced on a charged polymer in a polar binary mixture, 
(ii) In near-critical binary mixtures with salt 5 - 1 - and in 
polyelectrolyte a 52 i 53 ' 54 , ion-induced aggregates have been 
observed, where relevance should be the preferential sol- 
vation, (iii) We may replace hydrophilic counterions by 
hydrophobic ion o 39 i 45 such as tetrarphenylborate BPI14 . 
For hydrophilic and hydrophobic ion pairs, the preferen- 
tial solvation can be much stronger than for hydrophilic 
ion pairs^. Sadakane et al£^ added NaBPh.4 into a bi- 
nary mixture to obtain mesophases. (iv) Ionic surfac- 
tants can be adsorbed to DNA even if their bulk concen- 
tration is very smal l 24 ' 56 . The resultant complex can be 
solubilizcd in organic solvents such as ethanol 5 ^. Kuhn 
et al. predicted that this adsorption can occur as a first 
order phase transition 5 ^. 

In future work, we will treat two parallel rods in mix- 
ture solvent, between which there can arise attraction 
mediated by the composition fluctuations. 



Acknowledgments 

This work was supported by Grants-in-Aid for sci- 
entific research on Priority Area "Soft Matter Physics" 
and the Global COE program "The Next Generation of 
Physics, Spun from Universality and Emergence" of Ky- 
oto University from the Ministry of Education, Culture, 
Sports, Science and Technology of Japan. The authors 
thank K. Yoshikawa, K. Nishida, T. Sumi, Y. Masubuchi, 
and Y. Yamasaki for valuable discussions. 

Appendix A: Calculations for one-component 
solvent without salt 

We show exact results for one-component solvent with- 
out salt^i 3 ^. The function F\(r,a) in Eq.(2.16) depends 
on f = r/R logarithmically. Depending on whether 
a < a* or a > a* it behaves as 

Fl ( r f^ = [sinhCBlnr - tanlT 1 B)]~ 2 (a < a*) 
B z 

= [sin(.Blnf -tan^B)]" 2 (a > a*). (Al) 
The parameter B is the solution of the equation, 

1 ~ 1b ° = coth[BA/ + tantT 1 B] (a < a*) 
B 

= cot [BM + tan" 1 B] (a > a*) , ( A2) 

where we may set B > 0. At a = a* , we have B = and 
F t (r, a*) = (1 - lnf)~ 2 for any M > 0. At r = b it holds 

Fi(b,a) - (£b<t-1) 2 tB 2 , (A3) 
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where — is for a < a* and + is for a > a* in the right 
hand side. For M » 1 we find B = 1 — £b<t for -<;> 1 
and _B = tt/M for > 1 in terms of q in Eq.(2.20). 

At r — R we obtain the large M behavior F\(R,a) = 
1 - (1 - £ B a) 2 for £ B a < 1 and Fx(R,a) = 1 for £ B a > 
1. We may rewrite this result in terms of the osmotic 
pressure II in Eq.(2.24). In the limit M — * oo, it follows 
the Manning limiting law for the osmotic pressure, 



n = Tn p {\- £ B a/2) (£ B a < 1) 
S Tn p /2£ B a (£ B a > 1), 

= a/irR 2 . 



(A4) 



where n p = a/irR 2 . For £ B a > 1, II sat- 
urates at T/2tt£ b R 2 . Thus a fraction of 1 - 
(£bO')~ 1 of the counterions are apparently local- 
ized around the rod (the Manning-Oosawa counterion 
condensation)^>i£&2£&i2 

We calculate f2 in Eq.(2.13). Since the electric field is 
given by E = Te~ 1 [d(\nFi)/dr — 2/r], the electrostatic 
energy becomes 



F e = TLa + TL£ B - 1 F e (£ B <T), 



(A5) 



where J- e (s) with s = £ B a appears in Eq.(2.27). From 
Eqs.(Al) and (A2) some calculations yielcU^ 



T e {s) = (l + B 2 )M + hi 
= (1 - B 2 )M 



2s 



(a < a*) 



In 



1 



2s 



1 + B 2 



-B 2 

{a > a*). (A6) 



For large M we obtain approximate expressions 
"2(1 ~s) 2 



Ms) 



s 2 M 



M 



+ ln 
21n(s 



M~ 



(9<-l) 



2 - s 

-1 + Af- 1 ) (?>-l), (A7) 



where q = M{s — 1). In deriving the first line we have 
used Eqs.(2.21) and (A3). We introduce M" 1 on the 
right hand sides to avoid the logarithmic divergence at 
s = 1. At <j = <7* we have = M - 2 ln(l + M). 

Appendix B: Charged rod in one-component 
solvent with salt 

Here we examine the counterion density and the degree 
of ionization in one-component solvent with added salt, 
which is completely dissociated into cations and anions. 
We assume that the cations from the salt are of the same 
species as the counterions from the rod. For example, we 
suppose the chemical reactions: 

-COOX ^ -COO" + X+ (rod surface) 
XC1 -> X+ + Cl" (salt in bulk), 

where X=H or Na. The cation and anion densities are 
denoted by n\ and ri2, respectively. We treat the anion 
density at r — R as a control parameter and write it as 
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FIG. 13: Normalized density of cations ci(r) = £Bb 2 ni(r) 
(left) and anions C2(r) = £Bb 2 ri2(r) (right) vs r/b with salt 
at £b&o ~ 1.5 and A = 0.192. The salt density is given by 
c B = l B b 2 n B = 10" 1 , 10~ 2 and 10" 3 . 
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FIG. 14: Degree of ionization a as a functions of salt con- 
centration c B £ B b 2 n B with salt on a semi- logarithmic scale for 
A = 0.192 and 0.0192. Here £ B a Q = 0.5 < 1.0 (left) and 
£ B a = 1.5 > 1.0 (right). 



n B . In this case, the electric field created by the positive 
charges on the rod is screened for r > k^ 1 , where 



{^bub) 1 ' 2 



(Bl) 



is the Debye wave number far from the rod. For r 3> 
fc^ 1 , the system is homogeneous with Ui(r) = 7i2(r). If 
R 3> fcn j we obtain the results in the limit R — > oo. By 
setting $(oo) = 0, we may write the ion densities as ni = 
nse~ e */ T and ri2 = nse e */ T , where ni(oo) = 712(00) = 
n B . We rewrite the Poisson-Boltzmann equation as 



dr 2 



1 d\e$ , 2 , /e$ 



(B2) 



The boundary condition is given by Eq.(2.6) at r = b. 

We numerically solved the above Poisson-Boltzmann 
equation and the ionization equation Eq.(2.14) to obtain 
equilibrium n\(r), ni(r), and a. In Fig. 13, we plot the 
normalized ion densities, 



ci (r) = £ B b 2 n x {r), c 2 (r) = £ B b 2 n 2 (r), 



(B3) 



which tend to c B = £ B b 2 n B for r 3> kp 1 . We set 
£bVo — 1-5 and A = 0.192. Here c±(r) increases and 
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C2<V) decreases near the rod. In Fig. 14, we show a as 
a function of cb for £b&o = 0.5 and 1.5. We find that a 
decreases with increasing the salt concentration eg. We 
also find a — > 1 as cb — > for £b&o — 0.5 < 1, while a 
does not approach unity as cb — > for Ib^o = 1.5 > 1. 
This is consistent with the behavior of a in the limit 
M — > oo in the salt-free case discussed in Subsec.IIB. In 
the limit cb — > (with R — oo), a approaches unity for 
£b&o < 1-0 (though this limiting behavior is not seen for 
A = 0.0192 in the left panel), while it tends to a constant 
in the range (^bCo)" 1 < a < 1 for Ib&o > 1.0 as can be 
seen in the right panel. 

Thus, with increasing the salt density, a is reduced and 
the electrostatic interaction is screened. Experimentally, 
with addition of salt, highly expanded polyelectrolyte 
coils have been observed to shrink, eventually resulting in 
precipitation of the chains (phenomenon known as "salt- 
ing out" of polyelectrolytes^) . 

Appendix C: Water adsorption around hy- 
drophilic ions at small water content 

We here present a statistical theory of water adsorption 
to hydrophilic ions^. The system is in the cylindrical cell 
b < r < R and < z < L. Each solvation shell consists 
of v water molecules, where v = 1, • • • , S with S being 
the maximum number. The binding energy is Tw v . If 
«)„ > 1, the adsorption can be significant even for small 
bulk water composition 4>b (see Eq.(C8)). 

The number of ionized monomers is Ni = La. The 
numbers of ^-clusters composed of v water molecules 
are v Ni. The total number of the hydrated ionized 
monomers is then (3Ni with 



(Cl) 



The fractions (3 V are determined by minimization of the 
free energy of the form, 



^ = JVo0B(ln^B-l)+JVj(l-/3)In(l-)8) 



v - Wu), 



(C2) 



where iVo = V/vo with V = ir(R 2 — b 2 )L being the volume 
occupied by the solvent. The total number of the water 
molecules is fixed as 



N </> B +N I J2vj3 l , = Noct>B, 



where <p B is the volume fraction without adsorption. 
From dF a jdft v = under Eq.(C3) we obtain 

A, = (l-/3)^e w », (C4) 
= l-l/[l + 5>£e™"]. (C5) 

For Ni <C Nq or for R ^> b, we may set <fi° — <f> <C <j> even 
if (3 approaches unity. We then calculate the excess free 
energy AF S = F s - TN <f>%Qxi<t>% - 1) due to the water 
adsorption. Some calculations give 



AF S 
T 



N cj> B H^b/^D+N^b-^b 

-NjHi + Y.^^- 



ln(l 



"/?)] 
(C6) 



Here we set hx{<j> B / '(j>%) — 4>b I <fi B — 1 in the first line to 
obtain the second line for Nj <C No. 

The formation of solvation shells around the counteri- 
ons may be calculated in the same manner. Let Tw' v with 
v = 1, • • • , S' be the binding energy of ^-clusters of wa- 
ter molecules. For monovalent counterions, we find the 
free energy decrease in the same form as that in Eq.(C6) 
with w v being replaced by w' v . For a sufficiently small 
water density outside the solvation shells, we may use 
the results of one-component solvents in Sec. II if Ao is 
replaced by 



A = Ao-ln 



s' 



b B e w - 



(C7) 



The dissociation constant Ko in Eq.(2.16) is changed to 
Kq = ^oi exp(— Ao). and the parameter A in Eq.(2.23) 
is replaced by A — 2irlBb 2 Ko- 

Let the maximum of w v jv and w' v jv be w s . Significant 
ionization enhancement occurs for 



i>s > exp(-w s ), 



(C8) 



(C3) where the right hand side is small for w s ^> 1. 
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